arXiv: 1501.04418v2 [gr-qc] 19Jun2015 


LIGO-P1400252-V3 


Template-space metric for searches for gravitational waves 
from the inspiral, merger, and ringdown of binary hlack holes 

Chinmay Kalaghatgi,* Parameswaran Ajith,^ and K. G. Arun* 

‘Chennai Mathematical Institute, Plot HI, SIPCOT IT Park, Siruseri, 603103 Tamilnadu, India 
‘‘International Centre for Theoretical Sciences, Tata Institute of Fundamental Research, Bangalore 560012, India 

Searches for gravitational waves (GWs) from binary black holes using interferometric GW detectors require the 
construction of template banks for performing matched filtering while analyzing the data. Placement of templates 
over the parameter space of binaries, as well as coincidence tests of GW triggers from multiple detectors make use 
of the definition of a metric over the space of gravitational waveforms. Although recent searches have employed 
waveform templates coherently describing the inspiral, merger and ringdown (IMR) of the coalescence, the 
metric used in the template banks and coincidence tests was derived from post-Newtonian inspiral waveforms. In 
this paper, we compute (semianalytically) the template-space metric of the IMR waveform family IMRPhenomB 
over the parameter space of masses and the effective spin parameter. We also propose a coordinate system, which 
is a modified version of post-Newtonian chirp time coordinates, in which the metric is slowly varying over the 
parameter space. The match function semianalytically computed using the metric has excellent agreement with 
the “exact” match function computed numerically. We show that the metric is able to provide a reasonable 
approximation to the match function of other IMR waveform families, such that the effective-one-body model 
calibrated to numerical relativity (EOBNRv2). The availability of this metric can contribute to improving the 
sensitivity of searches for GWs from binary black holes in the advanced detector era. 


I. INTRODUCTION 

Second-generation laser interferometric gravitational-wave 
(GW) detectors such as Advanced LIGO [1] and Advanced 
Virgo [2] are expected to start their first data-taking runs by 
2015 and are expected to reach their design sensitivities in a 
few years (see [3] for a discussion of likely observing scenar¬ 
ios). In addition, an advanced interferometric detector, called 
KAGRA [4], is being constructed in Japan, and another one 
in India, called LIGO-India [5] is expected to be built in India. 
Coalescence of binary black holes (BBHs) are among the most 
promising sources for the first direct detection of GWs. With 
the designed sensitivities of these detectors, anticipated detec¬ 
tion rate of BBH coalescences is ~ 0.4 - 1000 per year [6]. 
Observation of GWs from BBHs is expected to make signifi¬ 
cant contributions to our understanding of fundamental physics, 
astrophysics and cosmology (see [7] for a review). 

The expected gravitational waveforms from BBHs can be 
accurately computed from General Relativity using appropriate 
approximation techniques or numerical methods (see, e.g. [8- 
10], for reviews). Different “template” waveforms correspond¬ 
ing to different parameters of the binary (such as the component 
masses and spin angular momenta of the black holes), are cross- 
correlated with the data looking for correlations that exceed 
certain threshold, indicating the presence of a GW signal. This 
technique is called matched filtering. Additional signal consis¬ 
tency tests and multi-detector coincidence tests are employed 
to further assess the true nature of the signal (see [11] for a 
detailed discussion). Several searches for BBHs have been 
performed in the past using data from the previous science 
runs of LIGO and Virgo [12, 13]. The non-detection by these 
initial instruments is consistent with our expectation of the 
astrophysical rates of BBH coalescences [6]. 

Coalescence of BBHs typically involve three stages: In the 
early inspiral stage, the radial velocity of the black holes is 
much smaller than their tangential velocity, which itself is 
much smaller than the velocity of light (y^ <s: c). Gravita¬ 

tional radiation reaction causes the binary orbit to continuously 


shrink and eventually the black holes move with relativistic 
velocities and ultimately merge with each other. In the final 
ring down stage, the merger remnant settles into a KetT black 
hole by radiating a spectrum of quasi-normal modes. Accurate 
analytical models of expected GW signals from the inspiral 
(see [14] for a review) and ring-down (see [15] for a review) 
stages are available for the last two decades. The analytical 
approximation methods cannot be applied to the merger stage, 
where the gravity is strong and highly non-linear. Due to this, 
the first searches for GWs from BBHs were performed em¬ 
ploying templates either describing the inspiral stage [16, 17] 
or the ring-down stage [18, 19]. In the case of BBHs consist¬ 
ing of intermediate-mass (mj, m 2 > IOOMq), only the merger 
and ring-down stages of the coalescence will be observable in 
ground-based detectors. Such signals appear as bursts of gravi¬ 
tational radiation. Several searches have been performed in the 
past that look for short-lived excess power in the data [20-22]. 

Recent breakthroughs in numerical relativity [23-25] have 
enabled us to accurately compute the expected gravitational 
waveforms from the merger stage. Large catalogs of numerical- 
relativity waveforms have become available now [26-28]. This 
has lead to the development of several analytical models coher¬ 
ently describing the GW signals from the inspiral, merger and 
ring-down (IMR) [29-44]. Recent LI GO-Virgo searches for 
GWs from BBHs made use of these IMR templates [12, 13], 
which significantly improved the sensitivity of these searches. 

Since the parameters of potential the GW signals buried in 
the detector data are not known a priori, the data has to be 
cross-correlated with a “bank” of template waveforms corre¬ 
sponding to different parameters. A discrete set of template 
parameters has to be chosen such that, for any signal, there is 
always a template that is sufficiently close to it. At the same 
time, in order to minimize the computational burden of the 
search, it is desirable to keep the number of templates in the 
bank to a minimum. In order to address this covering prob¬ 
lem in template placement, a geometrical method has been 
developed [45-47]. In this method, templates are placed in the 
parameter space such that inner product between neighboring 
templates is fixed to a predetermined value (called the minimal 
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match [46]), say 0.97. This ensures that the loss of signal-to- 
noise ratio (SNR) due to the mismatch between the signal and 
the closest template in the bank is acceptably small, say 3%. 

This geometric formalism introduces the notion of a metric 
in the space of GW signals, which allows us to place templates 
in the parameter space employing lattice-based methods [48- 
51] or stochastic placement methods [52-55]. In addition, the 
multi-detector coincidence test [56] employed in these searches 
also requires the knowledge of the metric. The metric can be 
computed from the template waveforms using the formalism 
introduced in [46]. Metrics in the space of post-Newtonian 
(PN) waveforms describing the inspiral part of the coalescence 
have been computed in the past [46, 55, 57, 58]. In the past 
searches for GWs from BBHs, even those employing IMR 
templates, the PN metric was employed in the construction 
of template banks and in the multi-detector coincidence tests. 
Although the degradation of the SNR due to this choice of 
the metric was not drastic, future searches will greatly benefit 
from the knowledge of the actual metric in the space of IMR 
waveforms. 

In this paper, we compute the metric in the space of the IMR 
waveform family IMRPhenomB [40]. We identify a coordinate 
system in which the metric is slowly varying over the parame¬ 
ter space of interest, which is desirable in the construction of 
template banks. We show that the inner product between the 
waveforms (known as the match) computed using the metric 
has excellent agreement with the exact numerical computation 
of the match. We also show that the match function computed 
using the IMRPhenomB metric agrees well with that computed 
numerically from other IMR waveform families, such as the 
EOBNRv2 [30]. This metric can be employed in the construc¬ 
tion of template banks as well as in multi-detector coincidence 
tests. We expect that this will significantly contribute to im¬ 
proving the sensitivity of searches of GWs from BBHs in the 
advanced detector era. 

The rest of the paper is organized as follows: Section II pro¬ 
vides a brief overview of the geometrical approach employed 
in the template placement, introducing the notion of the met¬ 
ric. Section III describes the calculation of the metric for the 
IMRPhenomB family and discusses the coordinates in terms of 
which the metric is more or less uniform. A discussion of the 
results, in particular the comparison of the match function com¬ 
puted using the metric with exact match function, is provided 
in Sec. IV, while Sec. V provides some concluding remarks 
and lists the future directions. 


II. GEOMETRICAL APPROACH TO TEMPLATE 
PLACEMENT 


Here we provide a brief overview of the metric formalism 
originally introduced in [46] for laying down waveform tem¬ 
plates in the parameter space of compact binaries. A set of 
intrinsic and extrinsic parameters A = (Aextrinsic,-^intrinsic) pa¬ 
rameterize a gravitational waveform /;(/; A). The intrinsic 
parameters are parameters that are intrinsic to the source, such 
as the masses and spins of the compact objects, while the 
extrinsic parameters are those which depend on the relative 
location of the source with respect to the detector (such as the 
time of arrival fg of the signal at the detector and the phase of 


the signal ipQ at a reference time fg *. 

The match between any two waveforms h{f. A) and h{f, A H- 
AA) is defined as: 


M(A, A + AA) = maxAA„.,„.,, {kf. A), kf, A + AA)), (2.1) 
where {a, b) denote the noise weighted inner product: 


{a,b) = 2 



a(m(f) + b{f)a*if) 
Shif) 


( 2 . 2 ) 


where Shif) is the one-sided power spectral density of the 
detector noise, /g is the low-freqnency cutoff of the detector, 
and a “hat” denotes a normalized waveform: h = h/\\h\\ where 

\\h\\ = v^. 

The match function has its maximum value (Aimax = 1) at 
A A = 0. Taylor-expanding the match function about A A = 0 
up to quadratic order gives: 

+ (2.3, 


The mismatch 1 - A4(A, A + AA) can be thought of as the 
proper distance between points A and A -H AA (in the signal 
manifold), and can be written as 

1 - M ^ gij AA'AAf (2.4) 

This introduces the notion of a metric in the parameter space, 
defined as 

1 / \ 

2{dAA‘dAAj) 


The metric over the intrinsic parameter space can be calcu¬ 
lated from the Fisher information matrix by projecting it on to 
the subspace orthogonal to the space of extrinsic parameters. 
The Fisher information matrix is defined as: 

r,7 = ^ {dkf;A),djkf;A)), (2.6) 

where 5, denotes partial derivative w.r.t. the parameter T,. The 
metric over the three-dimensional space of intrinsic parameters 
can be computed from this as 

g = ri-rTr 3 'r 2 , (2.7) 


where Fi is the Fisher matrix over the intrinsic parameters, Fg 
the same over the extrinsic parameters, and F 2 that describe 
the cross terms. Assuming three intrinsic parameters (say, 
total mass, mass ratio and one spin parameter as described in 
Sec III A), and two extrinsic parameters (reference time and 
phase), Fi, F 2 , F 3 are defined, respectively, as 


’ll 

ri 2 

ri3 

’21 

r 22 

r23 

’31 

r32 

F 33 


F41 F42 r 43 
Tsi r52 r53 


F 44 F 45 
r54 Tss 


( 2 . 8 ) 


* In the case of binaries not exhibiting spin precession, if we neglect the effect 
of higher harmonics (no-quadrupole modes), it can be shown that all other 
extrinsic parameters are degenerate with the parameters ?o and and can 
be absorbed into these. 
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Similarly for a two dimensional metric corresponding to, say, 
the two mass parameters, Fi, r 2 , r 3 read 


Til Fn 


r3i r32 


F33 F34 

r2i r22 


r4i r42 

’ 

F43 F44 


respectively. 


(2.9) 


III. COMPUTATION OF THE METRIC IN THE SPACE OF 
THE INSPIRAL, MERGER, RINGDOWN WAVEEORMS 

A. IMRPhenomB waveform 


The IMRPhenomB waveform models GW signals from the 
inspiral, merger and ringdown phases of a coalescing black 
hole binary with non-precessing spins [ 40 ]. The waveform is 
described by three intrinsic parameters: the total mass M = 
mi + m2, the symmetric mass ratio rj = mim2lM^, and a single 
effective spin parameter^ s {mixi + m2Xi)IM, where1,2 
are the dimensionless spins of the two black holes and mi,2 are 
the masses of the binary components. The waveform is written 
as h(f) = A(f) where the amplitude A[f) is defined as 


A{f) = 


' [1 + 2 il fi 


/o < / < /l 

/l < / < /2 


Wr£{f,f2,0-) / 2</</3 

(3.1) 

where C - ^377 and d is the effective distance to the 

source The frequencies /i ,/2 are the transition frequencies 
between the inspiral, merger and merger, ringdown stages. The 
amplitude is zero below the cutoff frequency /o and above 
the cutoff frequency The other quantities appearing in the 
above expression are, /' s ///i, v = (;rM/)^^^, -C(/,/ 2 ,o-) 
is a Lorentzian function with width cr centered around the 
frequency / 2 . The normalization constants and Wr make the 
amplitude smooth over the transition frequencies /i and / 2 . The 
phenomenologically introduced parameters ei = 1 .4547;^:' - 
1.8897 and 62 = -1.8153;t' + 1.6557 model the amplitude 
of the merger part, while a2 - -323/224 + 45 It// 168 and 
a-i - (27/8 - 1 177 / 6 ) are the 1.5PN accurate post-Newtonian 
corrections to the inspiral amplitude. 

The phase of the waveform is given by 


'P(/) = 2;r/fo + 00 


122, rjiP 




k=2 


(3.2) 


where fo and 0o are the time of arrival of the signal and the cor¬ 
responding phase by to and <po, respectively. The phenomeno¬ 
logically calibrated parameters ij/u describes the phase evolu¬ 
tion of the binary and are given in Table I of [40] . 


^ Effective distance is a combination of the luminosity distance, the antenna 
pattern functions of the detector and the inclination angle of the binary, 
which determines the observed amplitude of GW signals in the detector. 


B. Choice of coordinates 


It is desirable to have a coordinate system in which the 
metric is (at least nearly) constant over the parameter space. 
For example, the approximation for the mismatch given in 
Eq. (2.4) assumes that the metric is nearly constant over the 
two points in the parameter space for which the match is com¬ 
puted, and hence is slowly varying over the parameter space. 
Additionally, knowledge of a coordinate system in which the 
local density of templates (which is proportional s/\g\) over 
the parameter space is nearly constant is greatly helpful for 
stochastic placement methods, since the new random “pro¬ 
posals” of the templates can be drawn from a simple uniform 
distribution [55]. The conventional “physical” parameters of 
the binary, such as (M, ri,x} do not form a nice coordinate sys¬ 
tem for this purpose, since -^/i^ varies by at least 4 orders of 
magnitude over the parameter space of interest (IMq < mi, m 2 , 
mi + m 2 > IOMq, -0.95 < X 0.95). A set of coordinates 
termed chirp times were introduced by [45] in which the non¬ 
spinning PN metric is slowly varying. This was generalized to 
the case of non-precessing spins by [55]. The square root of 
the determinant of the PN metric computed in this coordinate 
system is found to have variation < 30 over the “low-mass” 
region in the parameter space (IMq < mj, m 2 < 2QMq) [55]. 

We have found that in the PN chirp time coordinate system, 
the square root of the determinant of the IMRPhenomB metric 
has a variation of ~ 100 over the parameter space of interest 
(IMq < mi, m 2 , mi + m 2 > IOMq, -0.95 < x ^ 0.95). Here 
we introduce a new coordinate system {^ 0 ,^ 2 ,^ 35 ], which is a 
modified version of the PN chirp time coordinates. 


128 7/(7rM/o)5/3 ’ 


4 T] (ttM fo)^!^ ’ 

7r(17022-9565.9;r) 

4(7rM/o)2/3 


(3.3) 

(3.4) 

(3.5) 


It can be seen that ^3 are nothing but the dimensionless 

chirp times introduced in [47]: 


^0 = 27r/oTo, ft = -2nfoT3, (3.6) 

where tq and T 3 are the familiar Newtonian and 1.5PN chirp 
times [45,59], and ff is the low-frequency cutoff of the detector 
noise. 

The PN chirp time coordinates are judiciously chosen to 
make the corresponding PN coefficients in the phasing formula 
(at OPN and 1.5PN order) linear in these coordinates. For the 
case of nonspinning IMRPhenomB waveforms, the OPN term 
in the phase is the same as that of the PN waveform. Also 
the 77 - independent coefficient (test-mass limit) in the 1.5PN 
term is the same as that of the PN waveform (see Table 1 of 
[40]). Hence the corresponding (dimensionless) PN chirp times 
(ft, ft) make the metric nearly constant over the two dimen¬ 
sional space of mass parameters. But, when spins are included, 
the 1 .5PN coefficient of the IMRPhenomB phase has terms up 
to quadratic order in 7/ and;^'. Further, for ease of computations, 
the coordinates have to be invertible and the transformation has 
to be bijective. With these aims, we write fty as a product of ft 
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FIG. 1 : Variation of the square root of the determinant of the IMRPhenomB metric over the parameter space. The plots show the contours of 
^(log Igl - log Iglmin)! where IgUin is the minimum value of |g| in each panel. In the top panel, the metric is computed in the [M, r],x] coordinate 
system, while in the bottom panel the metric is computed in the {^ 0 ,^ 3 , ^ 35 ) coordinate system [see Eq. (3.3) for definition]. It can be seen that 
the variation of the quantity over the parameter space in the bottom panels is significantly smaller than that in the top panels. 


and the terms linear in tj and rjx in the 1.5PN order coefficient 
of IMRPhenomB waveform phase. Admittedly, our procedure 
of finding the new coordinate system is somewhat ad-hoc and 
there may he a better coordinate system in which the variation 
of the metric is even smaller. However, we believe that this 
coordinate system is adequate for the purposes of template 
placement and coincidence tests. 

The physical parameters can then be written in terms of our 
new coordinates as 


C. Computation of the metric 

The codes to compute the metric are available as part of the 
LALSimulation package, which is part of the LSC Algorithms 
Library [60] - the core software package used for GW data 
analysis by the LIGO-Virgo collaborations. Here we provide 
some details of the numerical implementation of the compu¬ 
tation of the metric in LALSimulation. The IMRPhenomB 
waveform amplitude described in Sec. (Ill A) can be rewritten 
as 


M 

X 


5 6 

327r2/o^o 

I 25 I 

17022 ^35 

9565.9 9565.9 


' 25^2 ' 


1/3 


(3.7) 

(3.8) 

Above, 

(3.9) 


Ai(f) 


In this coordinate system, has a maximum variation of 
< 40 over the parameter space of interest. Figure 1 provides ^ 2 (/) 
a comparison of the ^J\g\ as computed in the two different 
coordinate systems. 


A{f) = 


'Ai(/), 

< A2(/), 

A3(f), 


/l < / < /2 
/2 < / < /3 . 


(3.10) 


X A(A) = Cau 

^ Kk{\) /<'^- 2 )/ 3 . 


k=0 


f(A)X(/,/2,cr); ^ = CwJ-^ 


- 7/6 


(3.11) 
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where A s [M, rj^x) i® the set of intrinsic parameters Simi¬ 
larly, the phase can be written as 


'?(/) = Yj 


<i-5)/3 


(3.12) 


k=Q 


where 


ifik - 

Vs = 


128?? 

128/7 


forkVS, 


(3.13) 


This allows us to rewrite the Fisher matrix dehned in Eq. (2.6) 


1 3 

Z [ 5,'P, A, 5,31') + (5,A«, djA,) ], (3.14) 


a=l 


where 


5,Ai(/) = Z 


ilk-l)l6 


k=Q 

2 


diA2{f) = (A)/' 


ik-DH 


k=0 


diA^if) = ^(\)di£(\J) + dia\)£{X,f), 


5,W) = j]5,-^,(A)/' 


<*-5)/3 


(3.15) 


k=0 


The Fisher matrix is computed by numerically integrating the 
expression Eq. (3.14) where the derivatives are computed an¬ 
alytically. Note that, in this way, the derivatives have to be 
evaluated only once for one computation of the metric, improv¬ 
ing the efficiency of the computation. 

The Fisher matrix computed in physical coordinate system 
A s {M, ri,x} can be transformed to the modified chirp time 
coordinate system 9 = {^o,^ 3 ,^ 3 s) in the following way: 

r'^j'^rj, (3.16) 


where J is the Jacobian matrix of coordinate transformation 
Jik = dAildOk- The Fisher matrix F' in the modified chirp 
time coordinate system is used to compute the template-space 
metric using Eq. (2.7). 


IV. RESULTS AND DISCUSSION 

In this section, we discuss the results for the two- 
dimensional metric of the nonspinning IMRPhenomB wave¬ 
form and the three-dimensional metric of the non-precessing 
IMRPhenomB waveform. Metric calculations have been done 
assuming aLIGOZeroDetHighPower noise power spectral 
density of advanced LIGO [61] employing a low-frequency 
cutoff/o = 20 Hz. 


3 Note that o-i = 0 in Eq. (3.11). 


A. Metric of the non-spinning IMRPhenomB waveforms 

The nonspinning IMRPhenomB waveforms, described by 
(M, Tj, to, (po), can be easily obtained by setting the effective 
spin parameter= 0 in the waveform discussed in Sec. Ill A. 
We then compute the corresponding 4x4 Fisher matrix using 
Eq. (2.6) and obtain the two dimensional metric by projecting 
the Fisher matrix orthogonal to to and (po using Eq. (2.7). The 
quantity has a large variation over the parameter space 
while using (M, rj) as the coordinates (see Fig. 1). Hence keep¬ 
ing the template placement problem in mind, we use (^o>^ 3 )> 
the PN chirp times coordinates, in the future calculations. 

Figure 2 compares the ellipses corresponding to constant 
matches (0.97 and 0.99) obtained semianalytically within the 
quadratic metric approximation (black ellipses) with the con¬ 
tours of the match function computed numerically (colored 
contours). Rows denote different total masses 20Mq, 50Mq 
and 100 Mq and the columns denote different symmetric mass 
ratios 0.25, 0.16, 0.08 and 0.02. It is clear from the hgure that 
the semianalytical and numerical contours agree very well for 
a match of 0.99 but they differ slightly for a smaller match 
of 0.97. This is likely due to the inaccuracy of the quadratic 
approximation to the match function [see Eq. (2.4)] for larger 
values of parameter differences. The inaccuracy of the metric 
approximation is the largest for binaries with large, comparable 
masses, likely due to the fact that they have a small number of 
GW cycles in the detector band. Despite this, the agreement 
between the semianalytical and numerical contours is excel¬ 
lent, suggesting that the nonspinning metric computed here 
can be used for construction of template banks and for multi¬ 
detector coincidence tests in searches employing IMRPhenomB 
templates. 

The IMRPhenomB waveform has been calibrated against 
numerical-relativity simulations with mass ratio < 4, and is not 
expected to be faithful towards GW signals from BBHs with 
large mass ratios [40] . However, several more recent models, 
in particular the EOBNRv2 model, have been constructed and 
are expected to be more faithful for large mass ratios as well. 
Having verified that our metric is able to produce an excel¬ 
lent approximation to the match function of the IMRPhenomB 
waveform, we now investigate the ability of the IMRPhenomB 
metric, and some variants of the PN metric, to model the match 
function of the EOBNRv2 waveform. 

The results are displayed in Fig. 3. It can be seen that the 
IMRPhenomB metric produces a reasonable approximation to 
the size of the 0.97 match contour at all points in the param¬ 
eter space; but the orientation of the ellipses start to deviate 
from that of the constant match contours for small mass ra¬ 
tios [t] —> 0). We also compare the numerical match contours 
of the EOBNRv2with two variants of the PN metric: 1) the 
metric of the restricted PN waveform with 3.5PN accurate 
phasing terminated at the last stable orbit frequency /lso of 
the Schwarzschild geometry, and 2) the metric of the restricted 
1.5PN waveform terminated at 1024 Hz. The former is the 
most accurate PN metric available now and the latter is the 
metric used for two of the previous searches (employing IMR 
templates) for GWs from high-mass BBHs on the LIGO-Virgo 
data [12, 13]. It is evident that the 3.5PN metric overestimates 
the size of the match contour significantly. The overestima¬ 
tion becomes severe as we go to higher mass systems. Using 
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T] = 0.25 


rj = 0.16 


n = 0.08 


T] = 0.02 



FIG. 2: Comparison of the match ellipses computed from the non-spinning IMRPhenomB metric (black ellipses) with contours of the “exact” 
match function of IMRPhenomB computed numerically (color contours). The rows correspond to different total masses (20Mq, 50Mq, IOOMq) 
and the columns correspond to different symmetric mass ratios (0.25,0.16,0.08,0.02). The solid (dashed) black ellipses correspond to a match 
of 0.99 (0.97) 


this metric in the template placement will cause severe under¬ 
coverage of the parameter space. Surprisingly, the metric of 
1.5PN waveform terminated at 1024 Hz provides a reasonable 
approximation to the EOBNRv2 match contours at low masses. 
However, in the high-mass regime (M > SOMq), the analytical 
match ellipses computed using this metric are significantly 
smaller than the exact match contours. This will cause signifi¬ 
cant over-coverage of the template bank when used in template 
placement. However, we note that it might be possible to 
improve the agreement of this metric (which, we believe is 
accidental) by tuning the termination frequency. 


B. Metric of the non-precessing spin IMRPhenomB waveforms 

The three-dimensional metric corresponding to the two 
mass parameters and the effective spin parameter can be com¬ 
puted starting from a 5 x 5 Fisher matrix corresponding to 
{^ 0 , > to, M, where represent the new coordinates in¬ 

troduced in Sec. IIIB. Projection of the Fisher matrix orthog¬ 
onal to (fo, 0o) will give the three dimensional IMRPhenomB 
metric for non-precessing binaries. Ellipses corresponding 
to the two-dimensional slices of the three-dimensional ellip¬ 


soid is shown in Fig. 4, which are compared against the two- 
dimensional slices of the contours of the three-dimensional 
match function computed numerically for IMRPhenomB wave¬ 
form. Each column specifies the total mass, symmetric mass 
ratio and effective spin parameter that are considered. The 
solid black ellipses represent the contours of match 0.99 and 
the dashed black ones correspond to a match of 0.97. We find 
excellent agreement between the ellipses obtained from metric 
and the numerical contours for a match of 0.99. The agree¬ 
ment is not so good for a match of 0.97 as it is for a match 
of 0.99. Such a disagreement we have already encountered in 
the case of two dimensional ellipses and is argued to be due to 
the break down of the (quadratic) metric approximation to the 
match function. However, in general, the agreement between 
the analytically computed match function from the metric and 
the numerical exact match contours is excellent. 


V. CONCLUSION AND FUTURE DIRECTIONS 

In this paper we computed the template-space metric of 
the inspiral, merger, ringdown gravitational waveform family 
IMRPhenomB. We proposed a coordinate system (which is 
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T] = 0.25 


rj = 0.16 


n = 0.08 


T] = 0.02 


20 Mo 





FIG. 3: Comparison of the match = 0.97 ellipses computed from the non-spinning IMRPhenomB metric (black ellipses) with contours of the 
numerically-computed “exact” match function of the non-spinning EOBNRv2 waveforms (color contours). Also shown are the match ellipses 
from the metric of 3.5PN waveforms truncated at LSO (red ellipses) and the metric of the 1.5PN waveforms truncated at 1024Hz (black dashed 
ellipses). The rows correspond to different total masses (20Mo, 50Mq, IOOMq) and the columns correspond to different symmetric mass ratios 
(0.25,0.16,0.08,0.02). 


an adaptation of the PN chirp time coordinates) in which the 
metric is slowly varying over the parameter space of interest. 
The semianalytical match function computed using the metric 
has excellent agreement with the exact match fnnction of the 
IMRPhenomB waveform computed numerically. In addition, 
we have also shown that the analytical match function com¬ 
puted using the IMRPhenomB metric agrees well with the exact 
match function of other IMR waveform families, such as the 
EOBNRv2 over the entire parameter space of interest. This 
will potentially allow us to use the IMRPhenomB metric in the 
searches for GWs from BBHs that employ more accurate IMR 
waveform families as well (the ones that are already available 
and those under development). There is ongoing work to em¬ 
ploy this metric in the construction of IMR template banks 
(primarily using stochastic placement methods) as well as in 
multi-detector coincidence tests [62]. We anticipate that this 
will contribute to improving the sensitivity of searches for GWs 
from BBHs in the advanced detector era. 

We note that there is an alternative proposal to use “exact co¬ 
incidence” criterion between multi-detector triggers, whereby 
the same template bank is used in multiple detectors and two 
triggers are considered coincident only if they are detected 


with exactly the same parameters (see, e.g., [63]). While this 
method has its own advantages, it is quite possible that a given 
signal is ascribed different parameters in two detectors due to 
the presence of noise, potentially missing a coincidence. Thus, 
designing a coincidence criterion based on the real metric of 
the templates is likely to be advantageous. However, only an 
apples-to-apples comparison can quantify the advantage or 
disadvantage of these methods, and we leave this as fnture 
work. 
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